Variable range cotunneling and conductivity of a granular metal 
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The Efros-Shklovskii law for the conductivity of granular metals is interpreted as a result of a vari- 
able range cotunneling process. The cotunneling between distant resonant grains is predominantly 
elastic at low T < T c , while it is inelastic (i.e., accompanied by creation of electron-hole pairs on a 
string of intermediate non-resonant grains) at T > T c . The corresponding E-S temperature Tes in 
the latter case is slightly (logarithmically) T-dependent. The magnetoresistance in the two cases is 
different: it may be relatively strong and negative at T <C T c , while at T > T c it is suppressed due 
to inelastic processes which destroy the interference. 
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Introduction. The low-temperature conductivity of 
most granular metals (both three-dimensional samples 
and thin films) exhibits a typical insulating behavior, 
characterized by the Efros-Shklovskii law 



(1) 



In the samples with low room-temperature conductivity 
this law is observed in the whole range of T (from room 
temperature down to liquid helium temperatures) L, 
3]. Two ingredients are known to be necessary 
for the existence of the behavior in usual doped 
semiconductors with localized impurity centers: (i) soft 
"Coulomb gap" in the electron density of states, and (ii) 
long-range electron tunnelling between distant centers of 
their localization. The original idea on the building of the 
Coulomb gap as presented in Ref . Q| was recently adapted 
to granular arrays in Ref. . It was argued that the prin- 
cipal source for this gap is random background charges 
and the physical mechanisms behind these charges were 
discussed in detail. The feature (ii) is quite natural for 
doped semiconductors as it follows from exponential de- 
cay of wave-functions of localized electrons. It is much 
less trivial for granular media where each grain is typi- 
cally connected by tunnel junctions to its nearest neigh- 
bors only, therefore the very origin of long-range tun- 
nelling needs some special explanation. 

In this Letter we demonstrate that the variable range 
hopping in a granular metal involves the so called co- 
tunneling process Q (either elastic or inelastic). The 
elastic cotunneling is effective, if the temperature is low 
enough (namely, for T < T c ~ C^ 1 (S ■ Ec) 1 ^ 2 , where 
8 is the characteristic level spacing in the grains, Ec 
is the characteristic charging energy, and C ~ 10 is a 
large logarithmic factor, see below). At T > T c the con- 
ductivity is dominated by inelastic cotunneling processes. 
We directly show this for the case of the granular metal 
with poor room-temperature conductivity (small inter- 
grain conductances g = (h/e 2 R) -C 1). We expect that 
the same is also true for samples with moderately good 
conductivity (cf. e.g. 0,0) with properly renormalized 
Tes, but this more delicate issue will be discussed in 



a separate publication. Experimentally, relative role of 
elastic vs inelastic cotunnelling processes could be de- 
tected by the presence of noticeable low-field magnetore- 
sistance in the hopping regime: while elastic cotunnelling 
is expected to lead to negative magnetoresistance like it 
was predicted for doped semiconductors 0, 0, 0, 0] , 
inelastic cotunnelling is intrinsically incoherent and the 
whole effect of magnetic field upon conductivity is local- 
ized within individual grains and can, therefore, be only 
observed in very high fields H > 10T . 

Co-tunnelling as a key mechanism of low-temperature 
charge transfer was proposed |2j and extensively dis- 
cussed (cf. e.g. 14] for the review) with regard to 
transport via quantum dots. Quantum dot situated 
between two bulk metal reservoirs is characterized by 
its charing energy Ec and dimensionless conductances 
9r,l- Semiclassical "orthodox theory" of Coulomb block- 
ade [15| predicts exponential suppression of conductance 
through quantum dot at temperatures T -C Ec, i.e. 
G or th oc min(<7L,5/j) exp(— Ec/T), due to low probabil- 
ity of creation of real state of the dot with extra electron 
charge. Co-tunnelling process, on the other hand, occurs 
in the next order of perturbation theory in (small) tun- 
nelling amplitudes i_R,L, but does not contain exponential 
suppression factor exp(— Ec/T) since the dot state with 
extra charge occurs as virtual state only. One should dis- 
tinguish two kinds of co-tunnelling processes, elastic and 
inelastic ones. Elastic process occurs when tunnelling in 
and out of the dot deals with the same intra-dot electron 
eigenstate a, thus it leaves the dot in exactly the same 
quantum state as it was before it. On the contrary, in- 
elastic co-tunnelling leaves behind it an excited electron- 
hole pair (since one electron tunnels in the dot and popu- 
lates some eigenstate h, whereas another electron tunnels 
out of the dot from another, p-th eigenstate). Elastic 
co-tunnelling contribution to the conductance scales as 
Gd oc gLgnS/Ec, whereas inelastic co-tunnelling contri- 
bution is Gind oc g L g R {T/E c ) 2 , cf. HQ- Thus, upon 
temperature decrease first the inelastic co-tunnelling pre- 
vails over classical conduction given by Gorth, and then 
at T* ~ \fEc5 it gives way to the elastic co-tunnelling. 
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Below we generalize the above ideas to the situation of 
variable-range hopping in granular arrays. 

In the spirit of the standard variable range hopping 
(VRH) theory, we consider a charge transfer between 
two distant grains i and j with anomalously small en- 
ergies e%,Ej <C Ec of the "charged ground state" (i.e. 
the states with an extra electron, or an extra hole; for 
explicit definition of the energy e and detailed discussion 
of the corresponding "density of ground states" see @|). 
Such a transfer between distant grains proceeds via a 
string of intermediate grains, where typically the energy 
of a state with an extra electron is high (e ~ Ec). As in 
the case of a standard single-particle tunneling via reso- 
nant impurity states, the entire process is realized as a 
coherent sequence of local hops between adjacent grains 
in the string. There is an important difference, however. 
For the single-particle problem, where the same electron 
has to tunnel sequentially through all the impurities in 
the string, (starting from the first one and ending with 
the last), the order of these local hops is fixed. In our 
case there are many electrons in each grain, all these 
electrons being ready to tunnel to an adjacent grain at 
any time, so that the sequence of local hops in such a co- 
tunneling process can be arbitrary. As a result, the states 
with a number of excited electrons and holes on differ- 
ent grains appear as intermediate virtual states of the 
cotunncling process, and the number of charged grains 
in these intermediate states can be larger than one. In 
the final state, however, all the charges of the interme- 
diate grains should be compensated, and the only extra 
charge is transferred between the two terminal grains of 
the string. Long-range hopping process involves many in- 
termediate grains; in general co-tunnelling through some 
of them will be of elastic type, whereas some other will 
be inelastic. We will see below that elastic cotunneling 
dominates the variable-range hopping at rather low tem- 
peratures T < T c , where T c is significantly lower than in 
the case of single intermediate quantum dot. 

Note also that in the case of not extremely low temper- 
atures T 3> S (when the spectrum of electrons in grains 
can be treated as a quasicontinuous one), one does not 
need to invoke phonons to ensure the energy conserva- 



tion: the energy can be taken from the fermionic ther- 
mostat via inelastic cotunneling process. 

General approach to variable-range cotunneling. The 
Hamiltonian of the system has the form H = Ho + H tun + 
He- Here the" single-grain" hamiltonian Hq = J^-Hq , 
and the inter-grain tunneling hamiltonian H tun = 
Srtj) ^tuni m the latter the summation runs over the 
pairs (ij) of neighboring grains, 

Mi) Vf n + a - V t n + a 

11 — / j c a i U' a . a u a i(n -"tun — / J L a i aj'J' ai a u 'ajU 

where the operator a^. a creates an electron in a single- 
particle orbital eigenstate ct^ with a spin projection a on 
a grain i. The Coulomb interaction hamiltonian 

Bc = lY,4 ii) (fii- *) fa-*), ( 2 ) 

ij 

hi = J2 a - a fl Q iff ffl a,a — n\ being the operator of excess 
number off electrons at the grain i, n$ corresponds to 
the minimum of Hq + He- Variables qi (not necessarily 
integer!) are the so called background charges (in the 
units of e) . We will treat them as independent continu- 
ous random variables —1/2 < q L < 1/2 with symmetric 
distribution P(q). Two different limits should be distin- 
guished: (i) the case of strong charge disorder, when the 
background charges are large, so that qi is distributed 
homogeneously in the interval —1/2 < qi < 1/2; and (ii) 
the case of weak disorder, when the charges are small, 
so that go = (q 2 ) 1 ^ 2 <C 1 and the probability to have 
q = 1/2, related to the "bare density of ground states" 
at the Fermi-level (cf. 0), is small: P 1/2 = P(l/2) < 1. 
While the former case seems to be most appropriate for 
a naturally disordered granular material, the latter one 
may be relevant for high quality artificial arrays of quan- 
tum dots. 

If the tunnel matrix elements t aiaj are small enough, 
the rate of cotunneling between distant grains i and 
j can be found in the high order perturbation theory in 

^tun : 
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Here we have assumed that the i —> j transition is dominated by tunneling along a unique "string" - a 
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chain of neighboring grains, denoted by numbers k — 
0, 1, . . . , N, N + 1, so that = i is the initial grain; 
N + 1 = j is the final grain; each pair k, k + 1 are in con- 
tact. The possibility for several relevant strings to exist, 
and the effect of interference of their contributions will 
be discussed in the last section of this Letter. The sum- 
mation runs over all possible partitions V of the string 
into two subsets: {k} = {m} U {I}; on the grains {m} 
the elastic cotunneling (via a state h m ,p m = a m ) oc- 
curs, while at grains {1} the inelastic cotunneling with a 
creation of an electron-hole pair with quantum numbers 
hi,pi takes place. The summation over all eigenstates a m 
and pairs of eigenstates hi ^ pi is assumed. The energies 
£pi , e hi are measured with respect to the Fermi level. The 
spin- variable <To corresponds to the state po; variables tr; 
correspond to ^-components of electron-hole pairs. All 
other spin-variables are not independent because of the 
spin-conservation by tunneling hamiltonian. The spin 
summation gives the factor 2 L ( V \ where L(V) is a num- 
ber of inelastic grains in the partition V . The interference 
cross-terms between the processes with different a m are 
neglected because of violent sign-fluctuations of t} lkPk l . 

The " time-orderings" T are all the possible orderings 
{fci, ^2, • ■ ■ , k r , , /ctv+i} of the set of grains k — 0, . . . ,N 
(note, that there is no final grain N + 1 in this set!). The 
contribution to the composite transition amplitude H^J 
corresponding to particular T has the structure 



fr(ij) (<T\ fV(fcjv+lfcjv + l + l) fi /yffcjvfejv +1) fi 



fi A-(fcifci 
' y £1 t\m 



where the many-particle Green-function Q = [Hq + 
Hc]^ 1 ■ The function Afc(T) = 1, if in the ordering T 
the grain k — 1 comes earlier, than the grain k, and 
Afc(T) = — f otherwise. The sign-factor (— 1) K ^ arises 
due to permutations of fermionic operators. The Fermi- 
functions /p(e) take into account the Fermi filling fac- 
tors. The (5-function ensures the energy conservation, 
and Ay = Sj — e, — is the difference of energies of 
the initial an final state. The factors 



N 



q(t) = n 



H c {n(r)} + J2(zh h 

r' = l 



(4) 



are products of energy denominators, appearing in Q. 
Here rife are numbers of excess electrons on fc-th grain 
after p local hops; they can be found from the following 
recursion formula: 

!n k (r - 1) - 1, if k = k r , 
»fc(r - 1) + 1, if k = k r + l, (5) 
nk{r — 1), otherwise, 

while Ufe(0) = is the equilibrium distribution. 

Inspecting the expression J2J we see that the character- 
istic value einei of the " inelastic energies" e p i , e^i is con- 
trolled by the combination of the 5-function in J3J and 



the product of the corresponding Fermi-functions. As a 
result, ejnei ~ A^/i, where L (which is T-dependent) 
is the number of inelastic cotunneling events in the total 
process. Actually T <C Cinei Ec, so that, in particular, 
the dependence of Q{T) on the inelastic energies can be 
neglected. On the other hand, the characteristic value 
of e e i of the " elastic energies" e am is limited only by the 
energy denominators Q(T), so that e e i ~ Ec- 

Thus, performing the integration over e p i,ehi and e Qm , 
we obtain for the effective Miller- Abrahams dimension- 
less conductance gij between two distant grains i and j: 



9ij oc e" 
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2N N 
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L=0 



(2L + 1)! 



N—L 



-F, 



NL 



(6) 



Here 5 = E e~ Xn ^ E ° vl ' s > is the mean level spacing [y is 
the electronic density of states per one spin projection 
at the Fermi level in a particular grain, v is the grains 
volume, and Eo is an arbitrary energy unit). The char- 
acteristic Coulomb energy Ec = E e ln ( Ekk / Eo ' > (normally 
Ec ~ e 2 /K e ffO, where a is the average diameter of grains, 
and K e ff is the effective dielectric permeability of the ma- 
terial, see 6]). Finally, he mean tunneling amplitude 
t = E Q e^ 21a ^ tkk +^ 2 / E o\ where \t kk+1 \ 2 is the "coarse 
grained" (i.e., averaged over an interval of energies near 
the Fermi level, large compared to the level spacing, but 
small, compared to any other relevant scale) value of 
\t ak _ a ' | 2 . The "averaging of the logarithm" rule ap- 
pearing in the above definitions of mean values, arises as 
a result of self-averaging of large (~ N) number of similar 
independent random factors with identical distributions. 
Note also the presence of the spin-factors 2 in the mul- 
tipliers, corresponding to inelastic processes and absence 
of such factors for the elastic processes. 

The local activation energy for the ij hop is the 
combination of £j, Ej and Ec , standard for the hopping 
conductivity theory (see |5| for the explicit definition). 

The weight-function Fnl = J2v 8l,lcp)C(V), where 
C(V) are numerical coefficients, depending only on the 

partition V and on the explicit form of the charging en- 

(fefe'\ 

ergy matrix E c : 



K(T)-K(T') 



T,T> 



Hdi m S XmiT) x mi T')0 [-A m (T)e Qm ] Q{T)Q{T') (7) 



e rn = e m /Ec and Q(T) = Q(T)E]j being the dimension- 
less variables. 

The explicit form of the weight-function Fn,l for gen- 

(kk' ) 

eral Ec can not be found. However, as it is argued 
below, the asymptotics of F, relevant for the purely 
elastic and purely inelastic limits, are F N ,o w if and 
_Fjv.jv ~ A-2 , correspondingly. The numerical constants 
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Ai and A 2 are not known. As a result, 

Ni- 



gtj cx exp 



HO 



I 8tt 2 £; c 
'e 2 i 2g |A,j| 2 ' 



clastic, 



(8) 



inelastic, 



where g = Gh/e 2 = %ir 2 {t/5) 2 <C 1 is the average dimen- 
sionless conductance of a contact between two adjacent 
grains. N ote, that this definition of g differs from that in 
Refs-HHEi]- A-PP^g standard Mott-Efros-Shklovskii 
arguments to the random network with conductances 
we obtain Eq. (fTJ with T ES = C(T)E C 



C\ In 



C(T) = 



c\ In 



8ir 2 E c 
A l9 S 
16ir 2 E 2 c 
e 2 A 2 gT 2 C 2 



T<T C , 
T> To, 



(9) 



where the crossover temperature T c ~ y/EcS/C, and 
ci ~ 1 is an unknown constant, depending on the statis- 
tical geometry of the granular material. S ince Eq cx a 
and 5 cx a -3 , we conclude that roughly Tes oc a -1 and 
T c cx a -2 . 

It should be noted that the above consideration is jus- 
tified and the VRH regime is actual only if the charac- 
teristic length of the hop is large: N ~ (Ec/CT) 1 / 2 > 1. 
For N < 1 the Nearest Neighbor Hopping regime, char- 
acterized by the Arrhenius law a cx exp(— Ea/T) with 
the activation energy Ea ~ Ec should be observed. The 
crossover temperature between the two regimes is con- 
trolled by the intergrain conductance g; the NNH is likely 
to be found in the samples with very low g. 

Model of local repulsion: Mott law for granular ar- 
ray. The expression can be explicitly evaluated for 
the model case of the short range Coulomb interaction 



( kk' ) 

Ec — EkSkk' ■ Then dimensionless local energies of 
the charged states are = (Ek/Ec) [1/2 =F 1k]i an d 

£j = Ek Q — \qi\) sign(<jj). In this case J7J can be writ- 
ten as a product of single particle Green-functions with 
energies, depending on the local charge. As a result 



c{v) = n 



i 



i 



-G/m 



n 



i 



i 



E, 



E 



so that, for N ^> 1, when the number of similar factors 
in I jl 0(1 is large and an effective self-averaging takes place, 



Fnl — C^A^ L A^ 



(10) 



A, = e-^*-^, A 2 = 4e ln * 2 - 2 MV4-« 2 ) ; (n) 

Cjq being binomial coefficients. The constant A\ does 
not show any dramatic dependence on the strength of the 



random potentials: A\ — e 2 « 7.4 for strong charge dis- 
order, and A\ — 4 for a weak one. The constant A 2 = e 2 
for a strong disorder, while for a weak one A\ ~ q 2 <C 1. 
The reason is the destructive interference between two 
possible processes of the pair production: in the "e — h- 
process" the electron is created first and the hole is cre- 
ated the second, while in the "h — e-process" the order 
is inverted. As a result, for go ^ 1 the crossover temper- 
ature strongly depends on q : T c (q ) ~ (q £)~ 1 \/EcS. 
The growth of T c is saturated at T™ ax - (£)-* (E^S) 1 / 4 
for q < (8/Ec) 1 / 4 , when the energies e p i,€hi of the pairs 
come into play. Thus, we conclude that for the case of 
weak charge disorder the inelastic cotunncling is sup- 
pressed and the crossover between elastic and inelastic 
cotunncling is shifted to higher temperatures. 

Unfortunately, the result ((10(1 can not be generalized 
for the case of nonlocal interaction Ec ■ It can be shown, 
however, that F/v,o ~ , and -Fjv,jv ~ A 2 , with certain 
rcnormalized constants A\ and A 2 . Roughly, the reason 
is as follows (details will be presented elsewhere) : simple 
exponential form of the F-function holds for any "effec- 
tively short-range" interaction (not necessarily strictly 
local one), while for effectively long-range one the func- 
tional form of F can be changed dramatically. The clue 
is that, despite the long-range character of the Coulomb 
potential, the interaction of effective degrees of freedom 
in our case is the short-range one. Indeed, actual charge 
configurations, relevant to our problem, are those, gener- 
ated by local electronic hops between neighboring grains. 
These hops create local dipoles, and the dipole-dipole in- 
teraction decays with distance r as r -3 . 

With the explicit formula ((10(1 at hand one can perform 
the summation in J^J and find 



Ingij = Nij < In 



AigS 
%n 2 Ec 



\£j - Sj 
2Nu 



2Ai 



AiEqS 



where the function <f(z) is implicitly defined by the rela- 
tions <p(z) — 2y — ln(I — y), y 3 = (1 — y)z 2 . The function 
y(z) = L/N (with asymptotics y{z) ~ z 2 ! 3 at z <C 1) has 
the meaning of a relative fraction of inelastic events. 

The average number of grains in a string, connect- 
ing two distant grains, is proportional to the distance 



between them: Ni. 



c 2 n g 



i/d 



'ij , where n g is the concen- 



tration of grains, and c 2 ~ 1 is a geometric constant, 
depending only on the statistics of grains packing. Es- 
timates, made for c 2 for several two-dimensional models 
show that c 2 ps 1. Thus, we have arrived at the d + 1- 
dimensional percolation problem in the and Si space. 
The density of sites Vd+i — TigPi/2/ Ec m this space is 
the density of marginal grains, whose ground states are 



T ' 
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almost degenerate. The connectivity criterion reads 



CfajEilrijEj) < f, where £(ri,£i\rj,£j) 



T 



hi 



A ig d 



2A ? _ 



AiEcS 



.12) 



As usual (see 0), one should find a value £ = £ c corre- 
sponding to the first appearance of an infinite cluster of 
grains, connected according to the criterion l|12fl . Then, 
with the exponential accuracy, the global conductivity of 
the system a oc exp(— £ c ). 

The arising percolational model differs, however, from 
the standard VRH one (see dEH) 

due to an additional 
dependence on ry and £j, ej appearing in the argument of 
the function ip. However, since the variation of this func- 
tion on the relevant scale of it's argument is Sep ~ 1 -C C, 
the corresponding relative variation of £(r^, Ei\rj, Sj) is 
small and can be treated by the standard perturbational 
method (see 0). As a result, we obtain the Mott law 



a ~ exp{-(T M /T) 



i/(d+i) 



(13) 



with Tm-, which is slightly temperature-dependent: 



2P 



1/2 



c 



2A 2 



aT/T l .)=r-,\\n[^^.)- X (T/T.. 



(15) 



The universal percolation constant n c f» 5.7 for d = 3 
and n c « 3.5 for d = 2. The function x(z) is related to 

ip(z) by 



(16) 



where the averaging over the " percolation hypersurface" 
in the space of dimensionless energy £ and dimensionlcss 
coordinate s has the following explicit meaning: 



(F) 



JdCdCJd S F(C,C,s)5(l-s- 



ICI+IC'I+IC-C'I 

2 



pcrc 



ICI+IC'I+IC-C'I 

2 



/ d(d(' J dsS (l - s 
The asymptotics of the function xi z ) are 

6iz 2/3 for z < 1 



2(lnz + 6 2 + .--) forz>l, 



(17) 



where &i,6 2 are universal constants, depending only on 
the space dimensionality d. In particular, 6 2 = —1/6 for 
d = 2 and b 2 = -1/2 for d = 3. 

Experimentally the Mott law (|13|) is likely to be ob- 
served in materials with weak charge disorder, where the 



density of states at e = is very small because of the fac- 
tor P\/2 <C 1, and the Coulomb gap is irrelevant, except 
for extremely low T range. In such materials, however, 
Tm is very large (cf. QD) and the crossover between 
the Mott and the Arrhenius laws should take place at 
relatively low T . 

Magnetoresistance. In the entire above consideration 
we have neglected a possible interference between contri- 
butions of different strings (if any), connecting the same 
pair ij. The main reason for this approximation is the 
strong sign-fluctuations of the matrix elements tf lkPk _ 1 - 
even in the coherent (purely elastic) limit the signs of 
contributions of different strings to the composite ampli- 
tude of i — * j transition fluctuate independently. The 
interference effects, although irrelevant to the zero-field 
effects, are sensitive to magnetic field, so that they may 
be the source for a low temperature effect of a negative 
orbital magnetoresistance, similar to the one discussed 
in [ToL ITU T12I H3I for the case of VRH in conventional 
disordered semiconductors. The key component of this 
effect is the interference between the contributions of dif- 
ferent spatial paths, leading from i to j. In our case this 
means the existence of several strings i — > j giving com- 
parable contributions to gij. For a fairly homogeneous 
material, where all gkk+i are of the same order of mag- 
nitude one can expect such different strings to exist al- 
ready for Nij <^ 1. For a strongly disordered material 
(with exponentially large fluctuations of gkk+i) there is 
typically only one leading string for Nij < iV m i n , and 
only for iVy > N min several strings act in parallel. The 
crossover length N m i n is a function of the magnitude of 
fluctuations D = (In g — In g) 2 , the explicit form of this 
function is model-dependent, and we will not discuss it 
in the present Letter. For us it is only important that 
N min (D) ~ 1 for D ~ 1, and N min (D) > 1 for D > 1. 

Apparently, the magnetoresistance is controlled by the 
typical area S(Ni oop ) of the interference loop (iVioop being 
the "length" of this loop). Loops with Ni oop < N min 
are extremely rare, and can be ignored. For N\ oop > 
Nmin the scaling law S(N\ oop ) ~ N^ op , with an unknown 
exponent u < 2 holds. The problem of statistics of loops 
is closely related to the well-studied theory of directed 
polymers in a random field |l7| . 

One of the essential ingredients of our cotunneling pro- 
cess is the presence of inelastic events, which certainly 
destroy the interference and suppress the magnetoresis- 
tance. Namely, the interference between two different 
strings A and B is possible only if the cotunneling at 
all grains of A and B, which are not common for them, 
is elastic (see Fig^l. Since the relative fraction of in- 
clastic cotunnelings y depends on the temperature, so 
does the length M c \ ~ 1/y of a typical stretch on a 
string, containing only "elastic grains". It is just M e \, 
not the entire length iV of the distant hop, that should 
play the role of the effective length 7Vi oop of the interfer- 
ence loop. Clearly, for T > T c one has M e \ ~ 1, while 



6 




FIG. 1: Four different strings contributing to the i —* j tran- 
sition in a particular realization of array. For a shown parti- 
tion ("inelastic grains" {1} are depicted as filled circles, elastic 
grains {m} - as open ones) only two strings (^4 and B) con- 
tribute to interference effects. 



M e i - (T c /T) 2 / 3 > 1 for T < T c . 

Thus, we can conclude that for temperatures T > 



T c N m f^ 2 (D) < T c the orbital magnetoresistance 
is strongly suppressed (since the typical elastic stretch 
is shorter than N min ), while at T <C T mag it can be rel- 
atively strong: the characteristic magnetic field H c , at 
which the conductivity would saturate at a(H 3> H c ) ~ 
V2a(H = 0) (cf. Ref.E3), is 



*o/5(M i)~ H^{T/T c ) 2u / 3 , 



(18) 



where = <&orig d is the field, corresponding to a flux 
$0 through an elementary triangle of neighboring grains. 
The dependence <j(H) at H < H c can be different in 
different r ang es of H: either Act oc H 2 (at the smallest 
fields, s ee lllL ll3j) or Act oc H (at the intermediate fields, 
see [ldLUl), or ^ a K H 1 / 2 (at relatively high fields, see 
0). 

Conclusions. In conclusion, we have developed a the- 
ory of variable-range hopping in granular arrays with 
poor intergrain coupling. Long-range hopping of elec- 
trons is provided by the multiple co-tunnelling "strings" 
which contains both elastic and inelastic processes within 
individual grains. In the presence of long-range Coulomb 
interaction, Efros-Shklovskii law for the temperature de- 
pendence of conductivity is derived in the asymptotic 
limits of purely inelastic or purely elastic co-tunnelling. 
Upon temperature decrease, relative contribution of elas- 
tic co-tunnelling increases; For the model case of local 
(screened) Coulomb interaction, general situation of par- 
tially elastic co-tunnelling was studied and crossover tem- 
perature T c was determined, cf. Eas. (|9I14|I . For real 
granular metals this crossover temperature happens to be 
rather low. In particular, for Al grains of size a ~ 20nm 
one estimates Ec ~ 500K and 5 ~ 0.05K, which leads 
(at g ~ 0.3, so that C ~ 12) to T c ~ 0.5K . For a ~ lOnm 
the same estimates give T c ~ 2K. Therefore the major 
part of experimental temperature range (from room to 
liquid helium temperatures) is dominated by inelastic co- 
tunnelling. This is the reason for magnetoresistance to 



be very weak in granular metals, contrary to disordered 
semiconductors. Observation of noticeable negative mag- 
netoresistance due to interference of different tunnelling 
"strings" might be possible with granular media made 
of small (< 10 nm) grains of non-superconductive met- 
als like copper, silver or gold, at temperatures below 
IK. In the inelastic regime (at T > T c ) the "constant" 
Tes, entering Eq£Q|, is itself T-dependent: it logarithmi- 
cally increases with the decrease of T. This dependence 
should lead to somewhat faster growth of the resistivity 
upon lowering T, than predicted by the standard Efros- 
Shklovskii law. 

After the present study was completed, we became 
aware of the preprint [lsj where purely elastic variable- 
range cotunnclling was proposed as the conduction mech- 
anism for granular metals; their results seem to agree 
with our ones as long as inelastic processes are neglected. 
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